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MODIFIED RAYLEIGH CONJECTURE FOR SCATTERING BY 
PERIODIC STRUCTURES 

ALEXANDER G. RAMM AND SEMION GUTMAN 



Abstract. This paper contains a self-contained brief presentation of the scat- 
tering theory for periodic structures. Its main result is a theorem (the Modified 
Raylcigh Conjecture, or MRC), which gives a rigorous foundation for a numer- 
ical method for solving the direct scattering problem for periodic structures. 
A numerical example illustrating the procedure is presented. 



1. Introduction 

For simplicity we consider a 2-D setting, but our arguments can be as easily 
applied to n-dimensional problems, n > 2. Let / : K — > M, f(x + L) = f(x) be an 
i-periodic Lipschitz continuous function, and let D be the domain 

D = {(x,y) : y>f(x), x EM.}. 

Without loss of generality we assume that / > 0. If it is not, one can choose the 
origin so that this assumption is satisfied, because M :— sup 0<x<L \f(x)\ < oo. 
Let x = (x, y) and u(x) be the total field satisfying 

(1.1) (A + fc 2 )u = 0, xeD, k= const >0 

(1.2) u = on S:=dD, 

(1.3) u = u + v, u : = e ikax , 

where the unit vector a = (cost?, — sin 8), < 8 < tt/2, and v(x) is the scattered 
field, whose asymptotic behavior as y — ► oo will be specified below, and 

(1.4) u(x + L,y) = vu(x,y), u x {x + L,y) = vu x {x, y) in D, v : = e lkL cos 6 . 

Conditions (|1.4|) are the qp (quasiperiodicity) conditions. To find the proper 
radiation condition for the scattered field v (x) consider the spectral problem 

(1.5) tp" + l 2 tp = 0, 0<x<L, 



(1.6) tp(L) = MO), V'{L) = V(0) 

arising from the separation of variables in (ll.Hl - l|1.4[l . This problem has a discrete 
spectrum, and its eigenfunctions form a basis in L 2 (0, L). One has 

if = Ae ilx + Be- ilx , A, B = const, 
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Ae ilL + Be~ ilL = v(A + B), HAe ilL - UBe~ aL = ilv(A - B). 

Thus 

p ilL _ u p—UL _ „ 

0. 



il{e ilL - v) -il(e- ilL - v) 

So, il{e aL ~v)[er aL -v) = 0. If 1 = 0, then tp = A + Bx, A + BL = uA, B = vB. 
Since v — e lkLcose ^ one has no eigenvalue I — unless fcisin^ = 27rm, to > is 
an integer. Let us assume that kLcosO ^ 27rm. Then 



AIL ikLcosO „ —ilL ikL cos 8 

e = e or e = e 



that is 

— — . or L =— /ccostH — 

L J L 



1+ = fccos6» + — -, or ^ = -fccosfl + — j = 0, ±1, ±2, . . . 



The corresponding eigenfunctions are x and e~ a i x . We will use the system 
e ll i x , which forms an orthogonal bais in L 2 (0, L). One has: 

/ e a tx e - a tx dx= e^ (3 ' m) dx = Q, j £ to. 
Jo Jo 

The normalized eigenfunctions are 

w (as) = — , j = 0, ±1, ±2, . . . 

These functions form an orthonormal basis of L 2 (0,L). Let us look for v(x) = 
v{x,y) of the form 

oo 

(1.7) v(x,y) = ^2 c 3 v j(y)fj(x), y>M, c 3 = const. 

j=—oo 

For y > M, equation implies 

(1.8) v'j + (k 2 - fyvj = 0. 
Let us assume that I 2 ^ k 2 for all j. Then 

(1.9) v j (y) = e i ^y, 

where, for finitely many j, the set of which is denoted by J, one has: 

(1.10) Hj = (k 2 -l 2 )^ 2 >0, if l 2 <k 2 ,jeJ, 
and 

(1.11) ^ = i{l 2 - k 2 ) 1 ' 2 , if l 2 >k 2 7 j^J. 

The radiation condition at infinity requires that the scattered field v(x,y) be 
representable in the form (|1.7|) with Vj(y) defined by (|1.9|) - (|1.11|) . 

The Periodic Scattering Problem consists of finding the solution to (|1 . 1|> — 
(|1.4|) satisfying the radiation condition l|1.7|l . (|1.9|l - (|l.ll|) . 

The existence and uniqueness for such a scattering problem is established in 
Section 2. Our presentation is essentially self-contained. In pQ the scattering by 
a periodic structure was considered earlier, and was based on a uniqueness theo- 
rem from [Jj. Our proofs differ from the proofs in pQ. There are many papers on 
scattering by periodic structures, of which we mention a few P, 0, @], 0, EH 
fTT)j.[TTj. [T2*] . [T3| . [T5] . |2Sj- The Rayleigh conjecture is discussed in several of the 
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above papers. It was shown (e.g. 0) that this conjecture is incorrect, in gen- 
eral. The modified Rayleigh conjecture is a theorem proved in |18| for scattering by 
bounded obstacles. A numerical method for solving obstacle scattering problems, 
based on the modified Rayleigh conjecture is developed in 8^. The main results of 
our paper are: the modified Rayleigh conjecture for periodic structures (Theorem 
4.4) and a rigorous numerical method for solving scattering problems by periodic 
structures, based on the modified Rayleigh conjecture (Section 4). The proof of the 
limiting absorption principle (LAP) and the rigorous and self-contained develop- 
ment of the plane wave scattering theory by periodic structures is also of interest 
for broad audience. This theory is based partly on the ideas developed in |21|. 
|22j . |19j . The proof of the key lemma 2.2 is based on a version of Ramm's identity 
(2.16). Numerical implementation of the method for solving scattering problems 
by periodic structures, based on the modified Rayleigh conjecture, is constructed 
using the approach developed in |SJ and in |2.'ij . Applications to inverse problems 
are discussed in [18] and |24*). 

2. Periodic Scattering Problem 

Existence and uniqueness of solutions of the Periodic Scattering Problem can be 
proved easily, if one establishes first the existence and uniqueness of the resolvent 
kernel G{x, y, £, 77, k) of the Dirichlet Laplacian in D: 

(2.1) (A + k 2 )G(x,y,tv,k) = -S(x-0S(y-v), G = on S, 

(2.2) G(x + L,y,£,ri,k) = vG(x, y, £, 77, fc), G(x,y,£ + L,rj,k) = VG(x,y,t; ) rj ) k), 
(2.3) 

G x (x + L,y,£,r),k) = vG x (x, y, £, 77, k), G x (x,y,£ + L,rj,k) = vG x (x,y,£,ri,k), 

and G satisfies the LAP, see (2.5) below. The overbar here and below stands for 
the complex conjugation. 

Indeed, if such a function G exists, then v can be found by the Green's formula 

(2.4) v(x,y) = - 7t (£, 77)67^(2;,?/, £,77,*;) ds, 

where N is the unit normal vector to S pointing into D. 

To prove the existence and uniqueness of G{x, y, £, 77, k) define 

io = -A 

to be the Laplacian on the set of C 2 (D) quasiperiodic functions vanishing on the 
boundary S, and vanishing near infinity. Let 

D L : = {(x,y) : < x < L, (x, y) e D}. 

Then Dl is a section of D, and £0 is a symmetric operator in L 2 {Dl). This oper- 
ator is nonnegative , and therefore [5] there exists its unique selfadjoint Friedrichs' 
extension, which will be denoted by I. 

Let Im(k 2 ) > 0. Then there exists a unique resolvent operator (£ — k 2 )^ 1 . Thus 
its kernel G(x, y, £, 77, k) also exists and it is unique. To establish the existence and 
uniqueness of the kernel for k > we are going to prove the following 
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Limiting Absorption Principle (LAP). Let k > 0, e > and assume that 
k 2 is not equal to A 2 . Then the limit 

(2.5) lim G(x, y, £,r),k + ie) = G(x, y, £, r], k), 

exists for all (x, y) £ D, x ^ y. The proof is based on the following two lemmas. 
Lemma 2.1. Let < e < 1, and a > 2. Then 

, 9R x f \G(x iy ,tv,k + te)\ 2 

(2-6) {1+e+v 2)«/2 d ^ * c ' 

where c = const > does not depend on e, and (x, y) is running on compact sets. 
Proof of Lemma \2.1\ It is sufficient to prove that the solution to the problem 

(2.7) (A + k 2 + ie)w e = F, in D L , w e £ L 2 (D L ), w e = on S L 

(2.8) w e (x + L,y) = i/w e (x,y), w ex (x + L,y) = i/w ex (x,y), 
satisfies the estimate 

(2.9) N*:=sup[ ^f' V ^ dxdy: = N 2 ( We )<c, 

0<e<ljD L (1+X^ + y^) a / 2 

where F S C^(Dl) is arbitrary, and c = const > is independent of e > 0. 

If l|2.9|l fails, then N en — > oo, e n — > 0. Define ip e : = w e /N t , where e : — e n . 
Then N(tp e ) = 1, ip e solves (j2~7l (with F replaced by F e := F/N e ), and satisfies 
(|2.8|) . From N(tp e ) — 1 it follows that ip e — ip as e — > 0, where — denotes the weak 
convergence in L 2 (Dl, 1/(1 + a; 2 + y 2 ) a/ ' 2 ) : = L 2 . By elliptic estimates, -0 e ^ ^ in 
h Io C ( D l), and therefore V> £ -> V in L 2 loc {D L ). This and (ETJl-lgSJl imply ip e -> ip 
in Hf oc (D^). Thus -0 solves the homogeneous (_F = 0) problem (|2.7|I - H2.8|I . If we 
prove that ip = 0, then we get a contradiction, which shows that l|2.9|l holds. The 
contradiction comes from the relationship = N(ip) = lim e ^ -^(V'e) = 1- O nc 
proves that 

(2.10) lim N(tp e ) = N(ip) 

e—>0 

as follows. If 

{x,y)ED R :={(x,y) : f{x)<y<R, < x < L}, 

where R > M is an arbitrary large fixed number, then lim^o N{ip t rju) = N^ipr/n), 
where 

(l, f(x)<y<R, 
[0, y>R. 

In the region D' R — {(x,y) : y > i?, < x < L}, one has \ip e (x,y)\ < c, (x,y) £ 
D' R . Thus 

sup N(iP e ( XL - %)) < O ( ) , < 7 < a - 2. 
0<e<l V-n-V 

The desired result (|2.10|l follows. 

To complete the proof let us show that the problem H2.7fl - I|2.8[l . with F — 0, and 

e = 0, has only the trivial solution w, provided that w is "outgoing" in the sense 

f L — 

Wj y — i^jUtj = o(l), as y — > oo, Wj : = / io^?j c?x. 

./o 
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One has 

(2.11) lim / (ww^ — w y w) ds = 0, 

Js R 

where Sr : — {(x, y) : y = R, < x < L}, ds = dx is the element of the arclength 
of Sr, and the overbar stands for the complex conjugate. 

Let us outline the steps of the further argument. 

Step 1: we prove that (|2.11(1 implies 

(2.12) w £ L 2 {D L ), \w\ + |Vw| < ce" 7lal , 7 = const > 0, 
if w is outgoing. 

Step 2: we prove that if w e L 2 (D L ) solves (2.7)-(2.8), with F = e = 0, then 
w = 0. Then we conclude that (2.9) (and (2.6)) holds, and, therefore, (2.5) holds. 
Let us prove (2.12). One has 



= / [w(A + k 2 )w - w(A + k 2 )w] dxdy 
Jd lr 



(2.13) = — I (wwn — wwn) ds + (wwn — ww^) ds 

Sl J Sr 

(wwn — wwn) ds, 

1 

where the Dirichlet condition l|2.7|l was used, and the integrals over the lines x = 
and x = L are cancelled due to the qp conditions l|2. 8fl : 



(-ww x + ww x ) dy+ {ww x - ww x ) dy 

x—0 J x—L 

(wwx — ww x ) dy — i vv{ww x — wwx) dy = 0. 

x=0 Jx=0 

Here we have used the relation vv = 1. Thus (|2.13l) implies 



(2.14) = / (ww y - ww y ) dx, VR > M. 

JS R 

If w is outgoing, then (|2.14() implies Wj(y) = for j 6 J, and \wj(y)\ < e~ 7 l y l , 7 = 
const > 0, so (|2.12|l holds. □ 

Lemma 2.2. Assume that w € L 2 (Dl), w solves \2.1^ with e = and F = 0, and 
w satisfies 1^2. tft) . Then w = 0. 



Proof of Lemma \2.2\ If w solves equation ()2.7|l with e = and F = 0, then w — 
w j(y) l Pj( x )- Since {<Pj(x)} is an orthonormal basis and w S L 2 (Dl), it follows 
that Wj(y) — for all j 6 J, and (|2.12l) holds. Let us use a version of Ramm's 
identity (JHI) p. 92), which is valid for any solution w of equation 11. 1|) which is 
outgoing in the sense that 

(2.15) w = '^^CjVj(y)ifij(x), Cj — const, j $ J. 



Note, that vj(y) = Vj(y) for j £ J . The identity is: 
(2.16) = (x 2 w,2W j ) J ■ + 1 \w t2 \ 
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where Wj :— dw/dxj, j = 1,2, x\ = x, X2 = y, over the repeated indices one 
sums up, |u>| 2 : = ww. The right-hand side of Ij2.16|l equals to 

- [X2(w,2jW t j — W.2jWj) + k 2 X 2 {w y2 W — W^w)] = 0, 

because wpw — w^w for outgoing w. 
One has 

(2.17) \w\ + \Ww\ < ce -7 ^! , 7 = const > 0, c = const > 0. 

Let R > maxf(x). Integrate (|2.1fi|l over Dlr '■ = {{x,y) : (x,y) € Dl, y < R} 
and use Green's formula to get: 



(k 2 \w\ 2 x 2 - \Vw\ 2 X2)N 2 



- _ v 2 

= — lira / [a^w^WjiVj + - — ! — ! — - — ^ ! — — — ] ds 

(2.18) JSlus r 

— lim / |^,2| 2 dx\dx2, 
where N is the normal pointing into Dlr, and we have used the relation 

(2.19) lim / \Vw\ 2 dxidx 2 = k 2 lim / |u>| 2 dx 1 dx 2 , 

which follows from the equation Aw + k 2 w — 0, boundary condition w = on S, 
quasiperiodicity of w, and from l|2.17[) . We have also used the relation w^wjNj = 
X2|Vw| 2 A^2, which follows from the condition u = on S. From l|2.18l) one gets: 

(2.20) lim / |w 2 | 2 dx x dx 2 = -- [ x 2 N 2 \Vw\ 2 ds. 
R ^°° Jd lr ' 2 J Sl 

Since f(x) is a graph, one has N2X2 > 0, and it follows from l|2.2U|l that w 2 = 0, 
soto= const, and const = because w\s = 0. Lemma 12.21 is proved. □ 

Remark 2.3. Condition of the type 

(2.21) N 2 x 2 > on S L 

was also used in [T§| . 

The proof of Lemma 2.2 is not valid if the Neumann boundary condition is 
imposed on S. 

3. Integral equations method 

In this Section we present another proof of the existence and uniqueness of the 
resolvent kernel G. We want to construct a scattering theory quite similar to the 
one for the exterior of a bounded obstacle ^7]. The first step is to construct an 
analog to the half-space Dirichlet Green's function. The function g = g(x, £, k) can 
be constructed analytically (x = (xi,X2),£, — {^1,^2))'- 

(3-1) S(x,f) = '%2<Pj(xi)(p j (£ 1 )g j (x2,£2,k), 

3 

9i :=gj{x2^2,k) = < . t 
to = (^y'c^ sin[^(6 + b)], fij = [k 2 - A 2 ] 1 / 2 , Vj (x 2 ) = e^ x \ 
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where 

ip'! + (k 2 - l])^ = 0, i>j{-b) = 0, Wfo, ipj] = 1, Aj = fccos(6>) + 

and W[v, ip] is the Wronskian. 

The function g is analytic with respect to k on the complex plain with cuts along 
the rays Xj — it, < r < oo,j — 0,±1,±2, in particular, in the region 5fc > 0, 
up to the real positive half-axis except for the set {A_y}j=o,±i,±2,...' 

Choose b > such that k 2 > is not an eigenvalue of the problem: 

(3.2) (A + fc 2 )V = 0, in£>_ 6 :={(x,y) : -b<y<f{x), < x < L}. 



ip\ y =~b = 0, ip N = on S, 

(3.3) 

ip(x + L,y) = uip(x,y), ip x (x + L, y) = vip x (x, y). 

One has 

(A + fc 2 ). 9 = -<5(x-0, x=(.t 1; x 2 ), ^=(6,6), 
xl ={(' T 7y) : — 6 < y < oo, < x < L}, 



(3.5) flr|„ — t = 0, 
and 

(3.6) (A + fc 2 )G=-<5(x-£),G = 0onS; 

G satisfies the qp condition and the radiation condition ( it is outgoing at infinity). 

Multiply H3.4fl by G, l|3.6l) by g, subtract from the second equation the first one, 
integrate over Dlr, and take R — > oo, to get 

(3.7) G = g+ (Gg N - G N g)ds = g - gfj, ds, n : = G n \ Sl . 

Js L Js l 

The qp condition allows one to cancel the integrals over the lateral boundary (x = 
and x = L), and the radiation condition allows one to have 

lim / (Gg N - G N g)ds = 0. 
R ^°°Js R 

Differentiate (|3.7() to get 

(3.8) n=-An+ 2^ on S L , Ay : = 2 jf ^fW) da. 

This is a Fredholm equation for /i in L 2 (Sl), if Sl is C fl,m , m > 0. The homoge- 
neous equation 1)3.8(1 has only the trivial solution: if [i + A/x = 0, then the function 
ip := f s gp* ds satisfies ipj^\s L — 0, where "0iv (V'iv) IS * ne normal derivative of 
ip from D-ij(Dl), and we use the known formula for the normal derivative of the 
single layer potential at the boundary. The ip satisfies also (|3.2f) and l|3.3|) , and, by 
the choice of b, one has ip = in £>_(,. Also ip — in D^, because (A + k 2 )ip — 
in Dl, ip\g L =0 (by the continuity of the single layer potential), ip satisfies the qp 
condition (because g satisfies it), and ip is outgoing (because g is). 

Since ip = in D-b and in Dl, one concludes that /i = ip~^ — ipjf, where ip^(ipjf) 
is the normal derivative of ip from £>_;, (£>£,), and we use the jump relation for the 
normal derivative of the single layer potential. 
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Thus, we have proved the existence and uniqueness of \x, and, therefore, of G, 
and got a representation formula 



(3.9) G = g- gfids. 

Js L 

This representation shows that the rate of decay of G as y — + oo is essentially the 
same as that of g. 

The G is analytic with respect to k on the complex plain with cuts along the 
rays Xj — it, < r < oo,j — 0,±1,±2, in particular, in the region Sfc > 0, up 
to the real positive half-axis except for the set {Aj} J= o,±i.±2,...- This follows from 
(|3.8(l . (|3.9(l . and the general result ^2], p. 57, [23j, concerning analyticity of the 
solution to a Fredholm equation with respect to a parameter. 

Suppose a bounded obstacle Do is placed inside Dl, u = on So = 9 Do, So is 
a Lipschitz boundary. If qp condition is imposed, then Green's function Go in the 
presence of the obstacle satisfies equations similar to (|3.9|) and (|3.8|l : 



(3.10) G {x,y) = G{x,y) - G(x, s)/j, (s, y) ds, fj, = G n, 

J So 

where N is the unit normal to So pointing into Dl, and 

,„ . . n dG f dG(s,a) , . 

(3.11) fi = -A p + 2^ on So, A O fi ■ = %J gN — fJ-o(cr) da. 

This is a Fredholm equation (with index zero). If k 2 is not an eigenvalue of the 
Neumann Laplacian in Do (=not exceptional), then equation l|3.11|l is uniqueley 
solvable and, by (|3.1U|I . Go exists and is unique for this k > 0. It is not known what 
are nontrivial sufficient conditions for k > to be not exceptional. The exceptional 
k form a discrete countable set on the positive semi-axis k > 0. If the Neumann 
boundary condition is imposed on Sl, then, even in the absence of the obstacle Do, 
it is not known if LAP holds, because the proof of Lemma 2.2 is not valid for the 
Neumann boundary condition on Sl. 

4. Modified Rayleigh Conjecture (MRC) 

Rayleigh conjectured [21] ("Rayleigh hypothesis") that the series ()1.7[) converges 
up to the boundary Sl- This conjecture is wrong (HS|) for some /(&). Since the 
Rayleigh hypothesis has been widely used for numerical solution of the scattering 
problem by physicists and engineers, and because these practitioners reported high 
instability of the numerical solution, and there are no error estimates, we propose a 
modification of the Rayleigh conjecture, which is a Theorem. This MRC (Modified 
Rayleigh Conjecture) can be used for a numerical solution of the scattering problem, 
and it gives an error estimate for this solution. Our arguments are very similar to 
the ones in [TH| . 

Rewrite the scattering problem 1|1.1|) - (|1.4|) as 

(4.1) (A + k 2 )v = in D, v = -u on S L , 

where v satisfies l|1.4fl. and v has representation l|1.7|). that is, v is "outgoing", it 
satisfies the radiation condition. Fix an arbitrarily small e > 0, and assume that 

(4.2) ||uo + V CjtfvjiyfojixiW < e, < x < L, y = f{x), 



where 



|j|<j(0 

U 2 (S L )- 



MRC FOR PERIODIC STRUCTURES 



9 



Lemma 4.1. For any e > 0, however small, and for any uo € L 2 (Sl), there exists 
j(e) and Cj(e) such that \4-ty holds. 

Proof. Lcmma l4~Tl follows from the completeness of the system {tpj (x)vj (f(x))}j=o,±i, 
in L 2 (Sl)- Let us prove this completeness. Assume that there is an h € L 2 (Sl), 
such that 

(4.3) / hipj{x)vj(f(x)) ds = 
for any j. From (|4 . 3|) one derives (cf. |17|. p. 162-163) 

(4.4) V>(x):=/ fts(x,0d£ = O, xeD_ 6 . 

Thus -0 = in L>i, and /i = ipt; ~ Vv ~ 0- Lemma 14. II is proved. □ 
Lemma 4.2. If 14-S\j holds, then 

|||u(x) - ^ Cj(e)?;j(y)^(a;)||| < ce, Va;,y S L> L , < x < L, y > f(x), 

lil<j'(<0 

where c — const > does not depend on e, x, y, and R; R > M is an arbitrary fixed 
number, and \\\w\\ \ = sup xeD \ DLR \w(x)\ + \\w\\ h i/2 {Dlr) . 

Proof. Let w : = v — Y^\j\<j( e ) c 3{ t ) v Ay) i PA x )- Then w solves equation jl.ljl . w 
satisfies I] 1.4(1 . w is outgoing, and ||HIl 2 (Sl) — e - One has (cf. I |2.4| l) 

(4.5) tu(x)=-/ wG N (x,0ds. 

Js L 

Thus i-e- IHI := \\w\\ L ^(s L ) < e, implies 

(4.6) \w(x)\ y=R < \\w\\ L 2 {SL) \\G N (x,Z)\\ L 2 iSL) < ce, c = const >0, 

where c is independent of e, and R > max/(x) is arbitrary. Now let us use the 
elliptic inequality 

(4-7) \\ w \\h™(d lr ) < c (11^11^-0.5(5^) + \\w\\ Hm - .s {SR) ) , 

where we have used the equation Aw + k 2 w = 0, and assumed that k 2 is not 
a Dirichlet eigenvalue of the Laplacian in Dlr, which can be done without loss 
of generality, because one can vary R. The integer to > is arbitrary if Sl is 
sufficiently smooth, and m < 1 if Sl is Lipschitz. Taking to = 0.5 and using (14.21) 
and 14.6fl one gets 

(4-8) \\w\\ h i/2 (Dlr) < ce. 

Thus, in a neighborhood of Sl, we have proved estimate (|4.8|) . and in a complement 
of this neighborhood in Dl we have proved estimate (I4.6|) . Lemma l4~2l is proved. □ 

Remark 4.3. In (|4.7|l there are no terms with boundary norms over the lateral 
boundary (lines x = and x = L) because of the quasiperiodicity condition. 

From Lemma f4. 21 the basic result, Theorem 14.41 follows immediately: 
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Theorem 4.4. MRC-Modified Rayleigh Conjecture. Fix e > 0, however 
small, and choose a positive integer p. Find 

(4.9) min \\u Q + Cj<Pj(n)vj(y)\\ ■ = m (p)- 

U\<p 

Let {cj(p)} be the minimizer of {4-^ - If m(p) < e, then 

(4.10) v(p) = J2 Cj{p)Vo{x)vj{y) 

\j\<p 

satisfies the inequality 

(4.11) \\\v - v(p)\\\ < ce, 

where c — const > does not depend on e. Ifm(p) > e, then there exists j = j(e) > 
p such that m(j(e)) < e. Denote Cj(j(e)) : = Cj(e) and v(j(e)) : = v e . Then 

(4.12) |||«-«e||| <ce. 



5. Numerical solution of the scattering problem 

According to the MRC method (Theorem 14. 4fl . if the restriction of the incident 
field — uq(x, y) to Sl is approximated as in l|4.9|l . then the series l|4.10|l approximates 
the scattered field in the entire region above the profile y = f(x). However, a 
numerical method that uses (|4.9|l does not produce satisfactory results as reported 
in |15| and elsewhere. Our own numerical experiments confirm this observation. A 
way to overcome this difficulty is to realize that the numerical approximation of 
the field —uq\s l can be carried out by using outgoing solutions described below. 

Let £ = (£i,6) € D- b , where b > 0, 

D-6 :={(&,&) : -6<6</(as), < < L}, 

and g(~K, £) be defined as in Section 3. Then g(x, £) is an outgoing solution satisfying 
Ag + k 2 g = in D L , according to (|3.4|) . 

To implement the MRC method numerically one proceeds as follows: 

(1) Choose the nodes jq, i = 1,2, ...,N on the profile Sl- These points are 
used to approximate L? norms on Sl- 

(2) Choose points C (1) ,C (2) ,-,^ (M) in D- b , M < N. 

(3) Form the vectors b = (uo(x;)), and a( m ) = (g(x 2 ,^ m ))), % = 1,2,..., A, m = 
1,2, ...,M . Let A be the N x M matrix containing vectors a( m ) as its 
columns. 

(4) Find the Singular Value Decomposition of A . Use a predetermined w m i n > 
to eliminate its small singular values. Use the decomposition to compute 

r mm =min{||b + Ac||, c e C M }, 

where 

t=i 

(5) Stopping criterion. Let e > 0. 
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(a) If r mm < e, then stop. Use the coefficients c = {ci, C2, cm} obtained 
in the above minimization step to compute the scattered field by 

M 

v{x,y) = c mg{x,y,t {m) )- 

m— 1 

(b) If r mm > e, then increase N, M by the order of 2, readjust the location 
of points G D_6 as needed, and repeat the procedure. 

We have conducted numerical experiments for four different profiles. In each 
case we used L — ir,k = 1.0 and three values for the angle 6. Table 1 shows 
the resulting residuals r rnin . Note that ||b|| = 1. Thus, in all the considered 
cases, the MRC method achieved 0.04% to 2% accuracy of the approximation. 
Other parameters used in the experiments were chosen as follows: A — 256, M = 
64, w m in — 10 -8 , b — 1.2. The value of b > 0, used in the definition of g, was 
chosen experimentally, but the dependency of r mm on b was slight. The Singular 
Value Decomposition (SVD) is used in Step 4 since the vectors aS m \ m = 1, 2, M 
may be nearly linearly dependent, which leads to an instability in the determination 
of the minimizer c. According to the SVD method this instability is eliminated by 
cutting off small singular values of the matrix A, see e.g. ^H] for details. The 
cut-off value w m i n > was chosen experimentally. We used the truncated series 
(|3.1|1 with \j\ < 120 to compute functions g{x,y,£/). A typical run time on a 333 
MHz PC was about 40s for each experiment. 

The following is a description of the profiles y = f(x), the nodes xj £ Sx, and 
the poles £ D-b used in the computation of g(xi, £( m )) in Step 3. For example, 
in profile I the x-coordinates of the A nodes Xj S Sl are uniformly distributed on 
the interval < x < L. The poles £< m ) e D_ b were chosen as follows: every fourth 
node Xi was moved by a fixed amount —0.1 parallel to the y axis, so it would be 
within the region £>_(,. The location of the poles was chosen experimentally to give 
the smallest value of the residual r mm . 

Profile I. f{x) = sin(2x) for < x < L, U = iL/N, x* = (U,f{U)), i 
1; 2, A, = (X4m,yim ~ 0.1), m = 1, 2, M. 

Profile II. f(x) = sin(0.2x) for < x < L, U = iL/N, jq = (U, f{ti)), i = 
1, 2, A, C (m) = (ita, Vim - 0.1), m = 1, 2, M. 

Profile III. f(x) = x for < x < L/2, f(x) = L-x for L/2 < x < L, U = 
iL/N, x, - (tj,/(ti))> ' l--' V. e (m) = (x 4m ,y 4m -0.1), m= 1,2,..., M. 

Profile IV. f[x) = x for < x < L, U = 2iL/N, x 4 = i = 

l,...,A/2, x, - (L,f(2(i-N/2)L/N)), i = N/2+1, ...,N,^ = (x 4m -0.03,y 4m - 
0.05), m = 1, 2, M . In this profile A/2 nodes x^ are uniformly distributed on its 
slant part, and A/2 nodes are uniformly distributed on its vertical portion x = L. 

The experiments show that the MRC method provides a competitive alternative 
to other methods for the computation of fields scattered from periodic structures. 
It is fast and inexpensive. The results depend on the number of the internal points 
£\ m ) anc j on their location. A similar MRC method for the computation of fields 
scattered by a bounded obstacle was presented in 
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